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Massive black holes are thought to reside at the centres of many galaxies™!, where 
[j they power quasars and active galactic nuclei. But most galaxies are quiescent, indi- 

cating that any central massive black hole present will be starved of fuel and therefore 
detectable only through its gravitational influence on the motions of the surrounding 
stars. M32 is a nearby, n quiescent elliptical galaxy in which the presence of a black 
hole has been suspectedutl; however, the limited resolution of the observational data 
and the restricted classes of models used to interpret this data have made it difficult 
to rule out alternative explanations, such as models with an anisotropic stellar ve- 
locity distribution and no dark mass or models with a central concentration of dark 
objects (for example, stellar remnants or brown dwarfs). Here we present space-based 
high-resolution optical spectra of M32, which show that the stellar velocities near the 
Q^ ■ centre of this galaxy exceed those inferred from previous ground-based observations. 

We use a range of general dynamical models to determine a central dark mass con- 
centration of (3.4 ± 1.6) x 10 6 solar masses, contained within a region only 0.3 pc across. 
This leaves a massive black hole as the most plausible explanation of the data, thereby 
strengthening the view that such black holes exist even in quiescent galaxies. 



We observed M32 with the Faint Object Spectrograph (FOS) on the Hubble Space Telescope 
(HST), on August 22, 1995. The COSTAR optics corrected the spherical aberration of the HST 
primary mirror. Spectra of the wavelength range 4572-6821 were taken through square apertures 
of 0.086" and 0.215", aligned on the nucleus and along the photometric major axis. The kinematical 
analysis was done by fitting each galaxy spectrum with the convolution of a template spectrum and 
a gaussian line-of-sight velocity profile. The details of the observations will be discussed by van der 
Marel, de Zeeuw & Rix ^manuscript in preparation) . The inferred kinematics are compared with the 
best ground-based dataP't3in Figure 1. The rotation curve derived from the HST spectra rises more 
steeply than that measured from the ground. The velocity dispersion from the HST aperture closest 
to the centre is 156±10km s _1 . The average of the four independent dispersion measurements inside 
the central 0.1" is 126 kins -1 , exceeding significantly the 85-95 kms -1 measured from the ground. 

Dynamical models are needed to infer the mass distribution from the kinematical data, and to 
test whether the observed velocities are consistent with a keplerian rise around a black hole. To this 
end we developed a new technique to construct fully general axisymmetric models that fit a given 
set of observations (Rix et al., van dec Marel et al., Cretton et al., manuscripts in preparation). It 
is based on Schwarzschild's approach", and can be-viewed as an axisymmetric generalization of the 
spherical technique of Richstone and collaborators^"". It consists of four steps: (i) an axisymmetric 
mass density is chosen that fits the M32 HST photometry" after projection; (ii) the gravitational 



potential is calculated, including a nuclear dark mass; (iii) a representative library of orbits is 
calculated; and (iv) a non-negative least-squares fit is performed to determine the combination of 
orbits that reproduces the assumed density and best fits the observational constraints. 

Previous studies of M32 have used either axisymmetric models with distribution functions 
depending only on the two classical integrals of motion (the energy E and angular momentum L z 
around the symmetry axis), f(E, L z ), or spherical models with fully general distribution functions. 
Models of the first type take into account the flattening of M32 (projected axial ratio ~ 0.73), but 
have a special dynamical structure. They successfully fit the ground-based data when a nuclear 
dark mass is includedlKI, but they leave open the possibility that more general models might fit the 
data without a dark mass. Models of the second type allow for general anisotropy, but ignore the 
flattening of M32. These spherical models were used only to interpret some of the older M32 datatl, 
and could not fit these without a nuclear dark massa However, we found that our less restricted 
axisymmetric models can easily fit these same data without a dark mass. So even though previous 
studies have shown that existing data are consistent with a nuclear dark mass in M32, it has not 
been shown that such a dark mass is required. 

To improve this, we used our general axisymmetric technique ton-determine which models can 
fit both the new HST data and the most recent ground-based data™^. The models take the ob- 
servational characteristics into account, and include velocity profile shape data (not available from 
the HST because of the limited spectral resolution of the FOS). Figure 2a shows a contour plot 
of the x 2 °f th e fit for edge-on models with a dark central point mass. The free parameters are 
the V-band mass-to-light ratio Ty of the stellar population, and the dark mass M 9 . The models 
labelled A-D are the best-fitting models for four fixed values of M,. Figure 3 compares the pre- 
dictions of these models with the observed major axis kinematics for each of the data sets. Models 
without a dark mass are ruled out; even the best-fitting model without a dark mass, Model A, 
fails to fit the nuclear HST velocity dispersions by > 40 kms -1 . Model C has M, = 3.4 x 10 6 M Q , 
and provides the overall lowest x 2 ■ However, we would not necessarily expect our algorithm to 
constrain a single, best-fit choice of the potential, since a three-dimensional distribution function 
might have the freedom to adjust to changes of the potential without affecting the fit to the data. 
Numerical experiments showed that the structure in the x 2 plot near the two minima could in fact 
be unreliable. We therefore do not claim that the local minimum labelled C is necessarily special, 
but conservatively estimate the allowed range for M, to be (3.4 ± 1.6) x 10 6 M . The uncertainty 
takes into account not only the observational errors in the data, but also the possible influence of 
small (< 2kms -1 in the projected kinematics) numerical errors in the models (due to, for example, 
gridding and discretization). M32 need not be edge-on, but models constructed for inclination 
i = 55° gave similar results and the same allowed mass range. A full exploration of-,the range of 
possible inclinations is computationally prohibitive, but our results and previous works suggest that 
the inferred M, depends only weakly on the assumed inclination. Allowing for possible triaxiality is 
not likely to change these conclusions; they are due mostly to the radial behaviour of the potential, 
and not to its shape. The best-fitting models are not f(E, L z ) models, but are similar in that they 
are tangentially anisotropic and require a similar dark mass. 

To constrain the size of the dark mass, models were built with the point-mass potential replaced 
by a Plummer potential with the same mass, but with scale length e. These were restricted to 
f(E,L z ) models for simplicity; constraints on e from the more general models should not differ 
by much. Figure 2b shows a contour plot of the \ 2 of the fit to the HST velocity dispersions, for 
edge-on models. The best fit has e = (a point mass); models with e ^ 0.06" are ruled out. 



The mass concentration thus has a half-mass radius r^ = 1.30e <J 0.08", implying a central 
density ^> 10 8 Mq pc -3 (1" = 3.39 pc) and a mass-to-light ratio ^ 20 r inside 0.1". It cannot be a 
cluster of ordinary stars; that would evolve rapidly by stellar collisions", and a change in mass-to- 
light ratio fromp£i|u2 at l" to ^> 20 at 0.1" would yield large broad-band colour gradients, ruled out 
by observations"". It must be dark, either a massive black hole or a cluster of dark objects (e.g., 
stellar remnants, brown dwarfs, planets). No theory predicts the formation of such clusters — except 
perhaps"" for a cluster of small black holes of mass m ~ 5-10 Mq, each — and most would not last 
for the age of the galaxy. A dark cluster is not an acceptable alternative to a massive black hole 
if it will collapse to a black hole in a short time, or expand or eject its mass through explosions or 
evolve in some other way that wouldptnake it observable (for then we would have to assume that it 
formed recently). Goodman and Lee— I considered possible dark clusters for M32, and argued that 
the clusters required r^ ^ 0.1" to be acceptable; these were consistent with ground-based data, 
which allowed" e to be as large as 0.4". We summarize and update the arguments here to show 
that all but the most implausible clusters are now ruled out. 



A cluster of stellar remnants of mass to ^ 1 Mq would undergo core collapse in a short time, 
5.5x 10 8 yr (MQ/m)(rh/0.08") 3 ' 2 for a Plummer model (and shorter for more concentrated models"). 
Neutron stars or small black holes would, form binaries and merge during the collapse by dissipating 
energy through gravitational radiation". The binaries would not be able to stop the collapse until 
they had grown large through multiple mergers, and even then they would-at most cause rapid os- 
cillations in the central density, during which they would continue to grow". One black hole would 
probably grow much faster than the others in a runaway manner". Its mass would be modest at first, 
perhaps 10 3 Mq, and might remain considerably smaller than 3 x 10 6 Mq. This model is not ruled 
out, but still requires a black hole (although embedded in a dark cluster) to fit the data. A cluster 
of white dwarfs would evolve similarly if the white dwarfs collapse to neutron stars upon merging; 
if they explode as supernovae instead the cluster would either lose mass or become observable. 

|A cluster of brown dwarfs or planets would have a long collapse time, but a short collision 
time". The outcome of a collision depends on the ratio of the collision and binding energies": 
brown dwarfs of mass m ^ 0.05 Mq would merge with little mass loss, leading again to the growth 
of a massive object; planets or small brown dwarfs would disintegrate, leading to a massive gas 
cloud that would be unlikely to remain dark. The constraint that the collision-time be longer than 
the age of the galaxy rules out clusters with m <^ 0.03 Mq. It might allow" a cluster of brow.ii 
dwarfs just below the hydrogen-burning limit (m ~ 0.08 Mq) if r^ ~ 0.08", but another argument" 
rules that out. The luminous stars of mass to* that pass through the cluster would get trapped by 
dynamical friction in a time that depends only on the cluster density and not on m (if m <C to*), 
5.5 x 10 8 yr (MQ/?7i i(! )(rh/0.08") 3 / 2 (the collapse time given above reduced by to/to*). The stars 
would sink to the centre, probably to merge into a massive object because binaries would be unable 
to eject stars from the deep potential well. A cluster with r^ ~ 0.08" would trap nearly 0.1% of the 
luminous stars, too many to be acceptable. This same argument limits the size of a bizarre dark 
cluster of elementary particles or tiny black holes to r^ <£ 0.01". 

The simplest and most plausible interpretation of the data is thus a single massive black hole. 
The X-ray emission from the centre of M32 could be from accretion onto the black hole, although its 
flux is low enough to be from-a. low-mass X-ray binary". The absence of a larger flux is explicable^: 
stars will be tidally disruptepp only once every 10 3 -10 4 yr, and the luminosity from continuous gas 
accretion can be suppressed". 



The arguments against a dark cluster in M32 rely mainly on the high density it would have 
(^ 10 8 M pc -3 ). These arguments are weaker for the active^alaxies M87 and NGC 4261, for 
which the dark matter density required by HST gas kinematicsE3El is only 10 5 (similar to the density 
of luminous material observed in some galactic nuclei and globular clusters). Stronger arguments 
can be made for three other black hole candidate galaxies requiring high densities: 10 7,3 for the 
quiescent galaxy NGC 3115, as derived from HST stellar kinematicsU: 10 9 ' 8 for our own Galaxy, as 
derived from near-infrared velocity measurements of individual-stars^; and io 9 - 6 " 12 - 6 for the active 
galaxy NGC 4258, as derived from water maser observationsS'Lj. Together with our M32 results 
this provides compelling evidence for massive black holes in both active and quiescent galaxies. 
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Fig. 1. — Stellar rotation velocities V (top panel) and velocity dispersions a (bottom panel) in M32, as 
function of major axis distance. The new HST data (aperture sizes-indicated in the figure) are compared 
with ground-based data obtained at the William Herschel Telescopet-I (WHT; resolution ~ 0.8"; connected 
by a line for illustration; error bars smaller than plot symbols) and the Canada-France Hawaii TelescopeP 
(CFHT; resolution ~ 0.5"). The HST data were obtained with the G570H grating on the FOS red detector, 
and were calibrated using arc lamp spectra obtained during Earth occultations. The pointing accuracy was 
limited by target acquisition uncertainties (~ 0.02") and thermal effects on the spacecraft and Fine Guidance 
Sensors (increasing to ~ 0.1" during the 14 hours of the observations). The actual aperture positiop-for each 
observation was calculated post hoc to ~ 0.01" accuracy from models based on HST photometry^-! for the 
target acquisition data, the observed count rate in each spectrum, and a verification image obtained at the 
end of the observations. The template spectrum used for the kinematical analysis was a mix of ground-based 
spectra of stars of different spectral types, chosen to best match the M32 spectrum. The kinematical results 
were corrected for line-spread-function differences between the galaxy and template spectra. The line spread 
function for each galaxy spectrum was calculated from models for the distribution of the galaxy light within 
the aperture, taking into account the grating broadening function and the finite size of the detector diodes. 
The relatively low velocity dispersion inferred for the HST aperture second-closest to the nucleus might have 
been affected by a systematic error, but none could be identified. 
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Fig. 2. — a, Contours of the \ 2 that measures the quality of the fit to the combined HST and ground-based 
kincmatical M32 data, for edge-on models with a fully general dynamical structure. The model parameters 
M. and Ty are the dark nuclear point mass and the V-band stellar mass-to-light ratio. The dots indicate 
models that were calculated. The predictions of the models labelled A-D are shown in Fig. 3. The contours 
were obtained through spline interpolation. The first three contours define the formal 68.3%, 95.4% and 
(heavy contours) 99.73% confidence regions. Subsequent contours are characterized by a factor two increase 
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The confidence values assume that there are only gaussian random errors in the data, 



and no numerical errors in the models. In reality there are numerical errors in the models, but they are 
known to be small (see the legend of Figure 3). Tests show that these errors can be responsible for the 
presence of the second (local) minimum in the x 2 contours (although this might also be real), but that they 
cannot make the predictions of models with M, < 1.8 x 10 6 M Q or M, > 5.0 x 10 6 M Q consistent with 
the data at the formal 99.73% confidence level. We therefore conclude that M. = (3.4 ± 1.6) x 10 6 M . b, 
Contours of the x 2 that measures the quality of the fit to the HST velocity dispersions, for edge-on f(E, L z ) 
models with an extended dark nuclear object with scale length e (and with Ty = 2.51 to best fit the data 
outside the central arcsec for this assumed dynamical structure). A point mass (e = 0) provides the best fit; 
models with e ;> 0.06" are ruled out. 
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Fig. 3. — Predicted rotation velocities V and velocity dispersions a of models A-D in Figure 2, compared 
with the major axis data from the HST, CFHT and WHT. Data points are plotted equidistantly along the 
abscissa. Corresponding major axis distances are given in the bottom panel. The WHT measurements are 
averages of data obtained at positive and negative radii. The n displayed quantities form only a subset of 
all the data that were fitted (for the ground-based observations"", velocity profile shape data and position 
angles other than the major axis are also available). Model C is the edge-on model with the lowest \ 2 , an d 
has a dark nuclear point mass M, — 3.4 x 10 6 M©. Models B and D are (approximately) the best-fitting 
models for M, = 1.9 and 5.4 x 10 6 M Q , respectively. Model A is the best fit without a dark mass. This model 
is ruled out. It is already marginally ruled out by the ground-based CFHT and WHT data, but only the 
HST velocity dispersion measurements make the case clear-cut. The accuracy of the model predictions was 
assessed by having them reproduce the known results for f(E, L z ) models"^. From these tests we conclude 
that numerical errors in our models are sufficiently small not to be of influence to the conclusions of our 
paper (they are < 2kms _1 in the projected kinematics and < 0.01 in the Gauss-Hermite momentg-J of the 
velocity profiles). 



